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Abstract 

We analyze the periodic tunneling of a Bose-Einstein condensate in a double- 
well potential which has an oscillating energy barrier. We show that the dy- 
namics of the Bose condensate critically depends on the frequency us of the 
oscillating energy barrier. In the regime of periodic macroscopic quantum 
tunneling (PMQT) with frequency usj, the population imbalance of the con- 
densate in the two wells can be enhanced under the condition of parametric 
resonance us = 2uij. Instead, in the regime of macroscopic quantum self- 
trapping (MQST), we find that MQST can be reduced or suppressed under 
the condition of parametric resonance between the frequency us of the energy 
barrier and the frequency usst of oscillation through the barrier of the very 
small fraction of particles which remain untrapped during MQST. 
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I. INTRODUCTION 



Macroscopic quantum tunneling with dilute Bose-Einstein condensates of alkali-metal 
atoms has been the subject of many theoretical [1-4] and experimental [5-7] studies. Here 
we consider the problem of trapped Bosons in a double-well potential and investigate the 
effect of a periodically varying barrier in the tunneling of the Bose condensate. We study the 
problem using the two- mode classical-like equations [1,2] and find that remarkable effects 
are related to parametric resonance [8] both in the periodic macroscopic quantum tunneling 
regime and in the macroscopic quantum self-trapping regime. 



II. TUNNELING IN A DOUBLE- WELL POTENTIAL 

In a dilute gas of N Bosons at zero temperature practically all particles are in the same 
single-particle state of the density matrix p(r, r'; t) [9]. This macroscopically occupied single- 
particle state is called Bose-Einstein condensate and its wavefunction ip(r, t) is well described 
by the Gross-Pitaevskii equation 
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where U(r) is the external potential, g = A7ch 2 a s /m is the inter-atomic strength with a s the 
s-wave scattering length, and the wavefunction is normalized to iV [10]. 

Let us consider a double- well external potential U(r). In the tunneling regime, the 
wavefunction of the Bose-condensate can be approximated in the following way 

V>(r,0 = jN^je^M*) + y/N^je^Mr) , (2) 

where Ni(t) is the number of particles in the first well and N 2 (t) is the number of particles 
in the second well, such as N = Ni(t) + N 2 (t) [1]. If the double- well is simmetric, it is not 
difficult to show that the time-dependent behavior of the condensate in the tunneling energy 
range can be described, with a suitable rescaling of time, by the two-mode equations 

c 



C = -Vl-C 2 sin0, = AC + - 7 === f cos0, (3) 



where ( = (Aq — N 2 )/N is the fractional population imbalance of the condensate in the two 
wells, = 0i — 02 is the relative phase (which can be initially zero), and A = 2E I / E T with 
Ei the inter-atomic energy and Et the tunneling energy, i.e. the kinetic+potential energy 
splitting between the ground state and the quasi-degenerate odd first excited state of the 
GP equation [1,2]. Note that E T = hu , where o;o is the oscillation frequency of the Bose 
condensate between the two wells when the inter-atomic interaction is zero (Ei = 0) [1-4]. 

The two-mode equations have been studied by many authors and they are the Hamilton 
equations of the following Hamiltonian 

H =^( 2 - V / W^cos0, (4) 

where ( is the conjugate momentum of the generalized coordinate 0. Actually, it is possible 
to introduce [1] an equivalent Hamiltonian, given by 

H' = \ V \ + aC 2 + bC , (5) 

where 

1 i A 2 

a = - v /l-AiJ , &= — , (6) 

with Hq = A£(0) 2 /2 — ^Jl — C(0) 2 cos 0(0). For the Hamiltonian H' the generalized coor- 
dinate is ( and its conjugate momentum is p^. The equations of motion derived by H' 
are 

C=P(, p f = -2«C-46C 3 . (7) 

The analysis of the two-mode equations (3) or (7) has shown that for A < 2 there is a 
periodic macroscopic quantum tunneling (PMQT) of the Bose-condensate with Josephson- 
like oscillations of frequency uj = u \/l + A [1,2]. Instead, for A > 2, there exists a critical 
C c = 2 V / A^T/A such that for < C(0) 

« ( c there is PMQT of condensate with frequency 
uj, but for ( c < ((0) < 1 there is macroscopic quantum self-trapping (MQST) of the 
condensate: even if the populations of the two wells are initially set in an asymmetric state 
(C(0) 7^ 0) they maintain, on the average, the original population imbalance with a very small 
periodic transfer of particles through the barrier with frequency ujst = Wq^JJ^Hq — 1) [1,2]. 
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III. AN OSCILLATING BARRIER IN THE DOUBLE- WELL POTENTIAL 



In this section we analyze the effect of a periodic oscillating energy barrier in the double- 
well potential. The presence of a oscillating barrier can be modelled by a time-dependent 
A. We choose the following form 

A(t) = A (l +e sin (cut)) , (8) 

where A is the static value, e is a small perturbation (e <C 1) and uj is the oscillation 
frequency of the barrier. 

We numerically solve the two-mode equations (7) with (8) by using a forth-order Runge- 
Kutta algorithm. First we consider the case of pulsed macroscopic quantum tunneling 
(PMQT). We set A = 0.8, C(0) = 0.5 and 0(0) = 0. In Figure 1 we plot the population 
imbalance (t, Q and its phase-space portrait (C? £>c) f° r some values of e and uj. For e = the 
population imbalance harmonically oscillates between the values —0.5 and 0.5 with frequency 
ujj = 1.27, but for e ^ the motion is not fully harmonic. In particular, Figure 1 shows 
that under the parametric resonance condition uj = 2ujj the dynamics of the population 
imbalance is modified: the amplitude of the oscillation is modulated and reaches the values 
— 1 and 1. 

Now we consider the case of macroscopic quantum self-trapping (MQST). We set A = 25, 
C(0) = 0.6 and 0(0) = 0. Figure 2 confirms that for e = the population imbalance 
does not change sign and it periodically oscillates between the values 0.47 and 0.6 with 
frequency ujst = 13.53. As shown in Figure 2, MQST is strongly affected by the parametric 
resonance between the frequency uj st of MQST oscillation and the frequency uj of oscillation 
of the energy barrier. Moreover, Figure 3 shows that at the parametric resonance condition 
uj = 2ujst and with a sufficiently large perturbation (e = 0.2) the system eventually escapes 
from the self-trapping configuration. Note that this remarkable effect has been confirmed 
[11] by using the nonpolynomial Schrodinger equation, an effective one-dimensional equation 
derived from the three-dimensional Gross-Pitaevskii equation [12]. 
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CONCLUSIONS 



We have studied the periodic tunneling and the quantum self-trapping of a Bose-Einstein 
condensate in a double-well potential with an oscillating energy barrier. We have used the 
two-mode classical-like equations to find that the periodic macroscopic quantum tunneling 
can be enhanced and the macroscopic quantum self-trapping can be reduced or suppressed 
by changing the frequency of the oscillating barrier. In the latter case, we have proved that 
the system escapes from the self-trapping configuration if the the frequency of oscillation 
of the double-well energy barrier and the frequency of MQST oscillations of the condensate 
satisfy the parametric resonance condition with a sufficiently large perturbation of the energy 
barrier. 
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FIG. 1. PMQT regime. Fractional population imbalance ((t) and its phase-space portrait 
(C,.Pc)- Initial conditions: C(0) = 0.5, 4>(0) = 0. A = A (l + esin(wi)) where A = 0.8. From 
top to bottom: (a) e = 0; (b) e = 0.25 and u = ojj; (c) e = 0.25 and oo = 2u>j; where ujj is the 
frequency of unperturbed PMQT oscillations. 
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FIG. 2. MQST regime. Fractional population imbalance ((t) and its phase-space portrait 
(C,Pc)- Initial conditions: ((0) = 0.6, (f)(0) = 0. A = A (l + esin(u;i)) where A = 25. From top 
to bottom: (a) e = 0; (b) e = 0.1 and lo = lust/2', (c) e = 0.1 and lo = 2ujst] where lost is the 
frequency of unperturbed MQST oscillations. 
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FIG. 3. MQST regime at the parametric resonance. Fractional population imbalance ((t) of 
the Bose condensate. Initial conditions: £(0) = 0.6, 0(0) = 0. A = Ao(l + esin (cot)) where Ao = 25 
and lo = 2ujst with lost the frequency of unperturbed MQST oscillations. From top to bottom 
and from left to right: (a) e = 0.01; (b) e = 0.05; (c) e = 0.1; (d) e = 0.2. 
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